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ABSTRACT 



The permanent beach erosion in Southern Monterey Bay is 
episodic, occurring infrequently when high tides coincide 
with stormy weather which allows wave action to erode the 
toe of the cliffs. Precise photogrammetric techniques are 
used to measure cliff recession from 1946 through 1984. This 
study shows maximun erosion occurs in the vicinity of Fort 
Ord (7.3 ft/yr) and decreases to the south. The analysis 
and errors associated with determining cliff recession using 
aerial photograinmetry are discussed in detail. 

A model is developed to predict cliff erosion based on 
the hypothesis that erosion only occurs when the water level 
due to combined tides, wave set-up and run-up exceeds the 
toe of the cliff elevation. The model combines predicted 
tidal elevations and wave heights. Shallow water wave 
heights at various locations are calculated by transforming 
deep-water directional wave spectra provided by the Fleet 
Numerical Oceanography Center. Refraction of the wave 
energy is responsible for the variability of erosion rates 
along the shore. The bathymetry of Monterey Bay is such that 
the refracted wave energy is greater in the Fort Ord area 
than to the south. The erosion model was calibrated using 
the spectral wave climatology and aerial photographs 
covering an 18 year period. The model qualitatively repli- 
cates the temporal variability of the measured recession 
rates and gives a reasonable prediction of the spatial vari- 
ation of the mean recession rates. 
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I. INTRODUCTION 



A. NEEDS FOR AND USE OF COASTAL EROSION DATA 

Most of the world's population is located in coastal 
areas, where the length of the world's coastline is about 
448,000 km. Residential and commercial buildings are 
frequently constructed close to beaches without regard to 
future physical changes of the coastline. In certain loca- 
tions the coastline is eroded by wave action. Often beach 
erosion leads to wave damage of coastal structures, which 
results in economic losses . 

These coastline changes, which are due primarily to 
wave action, illustrate the need for obtaining erosion data 
which can be used for rational coastal planning. Knowledge 
of shoreline . erosion rates is important for effective plan- 
ning and design on low sandy shores and on shores attacked 
directly by large storm waves. In such cases coastal 
erosion can be severe and permanent. Coastal development 
that is planned and constructed with appropriate considera- 
tion given to coastal erosion is the most economical. The 
need for obtaining erosion data is especially important for 
planning the growth of undeveloped areas of coastline. 
Erosion data may be obtained by studying coastal processes 
and collecting data over a long period of time using 
nautical charts, topographic maps, aerial photography, 
historical and repetitive field measurements. 

The southern shoreline of Monterey Bay, between the 
Salinas River and Monterey Harbor, is a sandy beach and 
backshore with dune fields. This shoreline represents one 
of the most scenic regions of Monterey county. Several 
structures have been built close to the shoreline without 
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apparent consideration of coastal erosion. For example, 
Stilwell Hall at Fort Ord was built 80 feet from the shore- 
line in 1946. It was necessary to construct a protective 
revetment in 1950; again in 1983, a temporary revetment 
had to be built after the shoreline had severely eroded. 
Similar problems have been confronted by the Ocean House 
apartments on Del Monte Beach where a temporary stone wall 
has been constructed to slow erosion. More importantly, 
there is considerable pressure to develop extensive sections 
of the coastline with new hotels, condominiums and apart- 
ments. Thus, for southern Monterey Bay , there is a crit- 
ical need for understanding erosion processes. 

B. COASTAL EROSION PROCESS IN SOUTHERN MONTEREY BAY 

Beach erosion is a process which occurs naturally in 
regions where the sea level is rising, as has been occur- 
ring during the past century along much of the California 
coast. Coastal cliffs are a common landform along large 
sections of the geologically young California coast. The 
cliffs are the result of active erosion, which occurs peri- 
odically during stormy periods . 

Monterey Bay, California's second largest, is semi- 
circular and opens to the west. The bay is 12 miles wide in 
an east-west direction .and 25 miles long in a north-south 
direction, (Figure 1.1). The outer limits of the bay floor 
correspond to the edge of the continental shelf. Monterey 
Bay cliffs consist of coarse sand which is highly suscep- 
tible to erosion from wave action and surface runoff. 

The sediment transfer and deposition budget along the 
shoreline depends on a source of sand, such as the Salinas 
River discharge and cliff erosion sand transport down the 
coast by waves and currents, and losses of sand into the 
deep ocean, onto land, and to sand mining (see Figure 1.2) . 



12 




Figure 1 . 1 



Monterey Bay 
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Figure 1.2 Sediment Transport Processes (Monterey Bay) 
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The difference between loss and gains determines whether 
there will be beach erosion or accretion. Studies conducted 
by Dorman (1968) and Arnal (1973) concluded that there is a 
net loss of sand, and therefore, a continuous erosion of 
Monterey Bay beaches. 

The waves of Monterey Bay consist of both Pacific Ocean 
swells with periods of approximately 12-20 seconds and 
locally generated wind waves. Swells are the most signifi- 
cant sources of wave energy into the Bay. Early spring and 
summer are typically characterized by swells from the north- 
west, while fall and winter bring swells from west and 
southwest (Figure 1.2). 

The distribution of sediments in Monterey Bay and its 
shoreline shape are in large part a response to wave 
climate. Sediment found near the shore, where wave effects 
at the bottom are stronger, tends to be coarser sand than 
that found offshore (Smith, 1983). Shorelines tend to 
orient themselves toward the waves that lead to the path of 
least resistance. The effect of the Monterey Bay bathymetry, 
anchored by rocky headlands at Point Pinos and Point Santa 
Cruz and perturbed by the Monterey Submarine Canyon, is to 
bend the wave- front into a shape approaching that of the 
shoreline shown in Figure 1.3 . 
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Wave Refraction Diagrams for Monterey Bay, (Wiegel, 1964 ) 




C. EARLIER STUDIES OF COASTAL EROSION IN SOUTHERN MONTEREY 



BAY 

Several coastal errosion studies have been conducted for 
Southern Monterey Bay; in general, they indicate a contin- 
uous erosion or recession of the coastal land. In most of 
these studies erosion is inferred by comparing shoreline 
changes measured from historical aerial photographs starting 
from the year 1939. Nautical charts from as far back as 
1856 and topographical maps are used to complement photo- 
graphs. Various reference points are selected for these 
studies to measure beach erosion, such as the top of the 
cliff, the toe of the coastal bluff, and the mean water 
line; selection of different reference points can lead to 
different measures of erosion rate. 

One of the more significant studies for our purpose was 
performed by Thompson (1981), who obtained consistent esti- 
mates .of erosion rates; he viewed aerial photographs with a 
mirror stereoscope to measure the recession of the toe of 
the cliff fronting the Phillips Petroleum Property. The 
average cliff recession rate between 1939 and 1978 was found 
to be 1.8 ft/yr (0.54 m/yr). 

Moffitt (1968) used aerial photographs to determine the 
erosion rate at the Monterey Sand Company; the reference 
point used was the waterline. Since the waterline depends on 
tide elevation, waves, seasonal beach profile, referencing 
the waterline results in a "noisy" measure of erosion rate. 
The average erosion rate found was 6.9 ft/yr (2.1 m/yr) 
between 1946 and 1967. Other studies have been performed in 
the area, and the resultig erosion rates from these studies 
are summarized in Figure 1.4 . 
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D. OBJECTIVES 



The objective of the thesis is to construct a model to 
predict coastal changes due to stormy weather, taking into 
account the major contributions to the coastal erosion 
process. It is hypothesised that erosion is episodic and 
occurs during simultaneous occurrences of high tides and 
storm waves; thus, it is proposed that cliff erosion due .to 
storms can be calculated, provided that local tide eleva- 
tion, wave height, and a function relating storm energy and 
erosion are known. 

The procedure followed is to obtain erosion rates using 
photogrammetry . Several of the previous erosion studies 
have utilized inprecise photo comparison techniques leading 
to erratic and questionable results. The precise photogram- 
metric techniques used here and associated errors are 
described in detail in Chapter II. The prediction of tides, 
the effect of waves and the various terms needed for stud- 
ying coastal processes are explained in Chapters III and IV; 
and the model which is used for prediction of coastal 
changes is described in Chapter V. 
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II. EROSION RATES MEASURED USING PHOTOGRAMMETRY 



A. METHODS FOR DETERMINING COASTAL EROSION 

Beach erosion rates can be determined by several 
methods. The best method to be used depends on -the avail- 
able data and the objectives of the proposed survey. 
Erosion rates can be measured by comparison of historical 
data obtained through old maps and nautical charts. The 
charts and maps indicate the location of the shoreline at 
the date when the field measurements were taken before 
compiling the map. The change of shoreline is determined by 
comparing the shoreline on the charts at different times. 
Annual recession rates for a specific coast are determined 
by dividing the measured change of the shoreline by the time 
interval. Considering that the shoreline recession in 
southern Monterey Bay is nominally 3 ft/yr (0.9 m/yr) and 
that the accuracy of standard 1:24000 scale USGS topo- 
graphic maps in horizontal positions is about 40 ft (0.5mm 
at scale), one would conclude that the measured recession 
rate using the above method is not accurate. Another 
problem is that the surveys for nautical charts are normally 
only made every fifty years in eighty percent of the areas 
concerned, and on five, ten, or twenty-five year cycles in 
the less stable areas that constitute the remaining twenty 
percent. Thus, only a small number of comparisons are 
available for most areas, with long time intervals between 
the possible charts ( Staff ord , 197 1 ) . 

Another method is to collect information for a specific 
beach through observations by local residents who have lived 
near the beach for long periods of time. Personal observa- 
tion can provide valuable information about changes, but it 
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is the least reliable method for obtaining accurate erosion 
data and is not generally useful for engineering planning 
purposes. This method is normally only useful in areas that 
have been developed for a considerable length of time. 

The most accurate method is to make repetitive profile 
surveys of the area of interest , utilizing permanent refer- 
ence points on the beach. However, this method is time 
consuming and expensive because the surveys must be repeated 
many times at relatively short time intervals. 

A less expensive alternative is to measure coastal 
erosion using aerial photographs. In the last four decades 
considerable progress has been made in pho to gramme t ry . 
Charting agencies presently use aerial photographs to 
produce charts. This procedure reduces field work and the 
compilation time. By utilizing photogrammetric methods the 
charts can be revised more frequently, and the data is more 
accurate. The photography of identical areas at short 
intervals makes it possible to observe changes in the 
shoreline. By drawing the shoreline imaged on aerial photo- 
graphs taken at different times, and comparing the drawn 
shorelines, coastal erosion rates may be estimated. 

Before utilizing one of the above methods the first 
question to be answered is: What location is indicative of 
permanent beach change? For example, does one measure 
changes at the mean waterline, at the high tide line, at the 
base of the backshore cliff, or at the top of this cliff? 
This question will be discussed further in Subsection 2.D.4, 
where it is described that the cliff tops are the preferred 
reference points for this study. 

The use of aerial photographs to determine shoreline 
changes has an advantage over other methods. For instance a 
permanent record of specific features, their location, and 
the contour of the shoreline are obtained in great detail at 
the time of the photography. All of the existing details of 
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the coast are recorded in the aerial photographs, while on 
topographic maps and nautical charts, details are limited by 
the purpose of the chart . For the purpose of measuring 
beach erosion, the aerial photograph permits a more compre- 
hensive analysis than maps or charts. Features such as high 
water levels, cliffs, waterlines, waves, reference build- 
ings, and roads are represented as a whole. In the last 
four decades governmental agencies such as the U.S. Coast 
and Geodetic Survey (USCGS), the National Ocean Service 
(NOS), the Defence Mapping Agency (DMA), and the U.S. 
Department of Agriculture (USDA) have been photographing the 
coasts more frequently than new maps or charts have been 
made. This allows the determination of coastal changes to be 
made more frequently and at relatively less cost than to 
produce a new map or chart (Stafford, 1971). 

The coast of Southern Monterey Bay has been photo- 
graphed at an average interval of about every five years 
since 1939. However, in the last decade the photographic 
interval has been shorter, about every three years, which 
allows a more precise determination of erosion rates and a 
clearer picture of coastal change. The process of taking 
the necessary aerial photographs is expensive; however, for 
this study, photographs were obtained from private and 
government agencies at relatively low cost ($150.00 per 
set). A disadvantage of the aerial photographic method is 
that aerial photographs are generally only available back to 
1939. 

Several errors may result from using aerial photography 
to measure coastal erosion. The most significant errors 
are : 

a. Scale variation due to altitude variation of the 
airplane . 

b. Scale variation due to ground elevation variation. 

c. Horizontal displacement of the image due to relief. 
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d. Scale variation due to tilt of the camera. 

A discussion of each of these errors is presented below. 

B. ERROR CONSIDERATIONS IN UTILIZING AERIAL PHOTOGRAPHS 

A map is a reproduction at reduced scale of an ortho- 
graphic projection of the terrain onto a reference datum 
plane. Each point on the terrain is seen to be projected 
normal to the reference plane. A distance measured between 
two points on the map when multiplied by the scale of the 
map will be equal to the corresponding horizontal distance 
measured in the field. 

A photo is a perspective or central projection. It is 
one in which all points are projected onto the reference 
plane through one point called the perspective center. A 
photo is only orthographic in the ideal case when the film 
and ground lie in parallel planes; therefore direct measure- 
ments on single photo may cause error (Moffitt, 1980). 
Because the study of coastal erosion requires very accurate 
shoreline determination, photogrammetric errors must be 
considered. The most important error inherent in photo- 
graphs is the scale variation. 

1. Scale Variation Due to Change in Ground Elevation 

The exact scale of the photographs must be known in 
order to obtain accurate ground distances from aerial photo- 
graphs. If ground distances are found by measuring photo 
distances on two different photographs and utilizing a 
nominal scale for the photographs , an apparent difference 
in ground distances can result where there is no actual 
difference. The scale, S, of a vertical aerial photograph, 
for a particular point is given by Equation 2.1 (Moffitt, 
1980) . 

S = f / (H-h) (2.1) 
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where f is the focal length of the aerial camera's lens, H 
is the flight height above datum, and h is the ground eleva- 
tion of the point. 




The flight height H is ordinarily controllable by the 
pilot to within 1 percent ( Slama , 1980 ) . Therefore, as can 
be seen from Figure 2.1a, if the photograph is vertical, 
and if there is little variation in ground elevation, there 
is only a small variation in scales of the photographs 
obtained in any given flight. Thus in Figure 2.1 b the 
scales for the points A and B are given by the equations 

S 4 = f/(H-h 4 ) (2.2) 

S B = f/(H-h 8 ) (2.3) 

where h 4 and hg are the ground elevations of the points A 
and B respectively. Specifications for aerial photographs 
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frequently require that if variations in scale exceed 5 
percent, the photographs be rejected which was a criteria 
used in this study. Scale variation can be controlled some- 
what by selecting appropriate flight/focal length combina- 
tion . 

The scale of an aerial photograph can also be 
obtained by comparing known distances utilizing the equa- 
tion 



S = ab/AB (2.4) 

where ab is the measured film distance between images of 
two points on the ground which are at the same elevation or 
where the difference in elevation is small, and AB is the 
distance measured on the ground (see Figure 2.1). 

A topographic map may be used to compute photo- 
graphic scales. This is not recommended since undefined map 
errors are then propagated into all the work done with the 
photographs. The ground distances between two points can be 
computed by measuring the map distance between points which 
are well defined on photographs, and multiplying by the map 
scale. In turn, the photo scale can be determined by meas- 
uring the photo distance between the same points, and using 
Equation 2.4 to calculate the average photo scale for one 
photograph; several scales for the same photograph must be 
computed for different ground points. A large-scale (rela- 
tive to the photo scale) map must be used to avoid measure- 
ment errors. The accuracy of this method will depend on the 
accuracy of the topographic map. 

A very accurate method for determining the scale of 
photographs is to utilize ground control points. The accu- 
racy of this method requires the existence of ground control 
points in the study area, good recovery of the existing 
ground control points, and images of ground control points 
that can be well identified on the photographs. 
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Ground control data are obtained using accurate 
surveying techniques. Ground control point identification 
is a most important step in photogrammetric mapping. Once 
the ground control data have been obtained, the photographic 
scale can be determined as follows: 

a. The photo distances between the control points 
are determined by measuring the photo coordinates of the 
images of the ground control points. 

b. Photographic scales are calculated utilizing 
Equation 2.4 each time for each pair of control points. 

c. The average photographic scale is taken as the 
average of these scales. By emloying an average photo- 
graphic scale, errors caused by tilt may be reduced. 

If the scale of a photograph is unknown, it can be 
determined by making use of other photographs with known 
scales . This method requires the existence of photographs 
with known scales which cover the same area as the photo- 
graph with an unknown scale and images of points well iden- 
tified in both photographs. 

The procedure for scale determination includes the 
following steps: 

a. The distance between two images is measured on 
the photo with the known scale. 

b. The ground distance AB between the two points 
is calculated using Equation 2.4 

c. The distance between the same images is measured 
on the photo with the . unknown scale. 

d. The scale of the photograph having the unknown 
scale is calculated from the Equation 2.4 

e. Steps a. thru d. are repeated several times for 
different image points to calculate the average scale of the 
photograph and improve accuracy. 

The accuracy of the computed photographic scale 
following the above procedure is not as high as for the 
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other methods mentioned previously because accuracy depends 
on the known scale of the photograph and elimination of 
other displacements . 

2 . Errors Due to Relief - Displacement 

As mentioned above, errors in measurements can be 
caused by elevation or relief differences. Because the 
ground is not usually flat, a photograph differs from a map. 
It is obvious from Equation 2.1 that the scales of the 
photographic images vary in accordance with the elevations 
of the corresponding objects on the ground, while the scale 
of a map is the same over the covered ground. Several 
methods for determining the average scale of a photograph 
were presented in the previous section. For a given 
method, the errors in the scale variation over the photo- 
graph can be reduced. 

Ground relief also affects the displacement of the 
images on a photograph. This displacement is known as 
relief displacement and is defined as the displacement of 
images radially inward or outward with respect to the photo- 
graphic nadir, depending on whether the ground objects are 
below or above the datum (Slama, 1980). 

The nature of relief displacement is demonstrated in 
Figure 2.2 . The ground points A and B are above the datum 
and their photographic images a and b have been displaced 
radially outward from the center of the photograph. If the 
ground points A and B were at the datum elevation, their 
images would appear at a' and b'. Relief displacement is 
zero at the center (o) of the photograph. By using the 
simple geometry of Figure 2.3, relief displacement can be 
determined from the following equations: 



d = r-r' 



rh/H 



(2.5) 



d = r'h/(H-h) 



( 2 . 6 ) 



where d is the relief displacement of the point, r is 
the radial distance on the photograph from the center to the 
image of the top of the object, r’ is the corresponding 
distance from the base of the object, H is the flight alti- 
tude, and h is the elevation of the object. 




Figure 2.2 Relief Displacement on a Vertical Photograph 

Equations 2.5 and 2.6, indicate that the magnitude 
of the displacement of an image is directly proportional to 
the elevation, h, of the object above or below the datum and 
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Figure 2.3 Relief Displacement of a Point 



to the radial distance r or r’ from the center of the 
image of the top or the base of the object respectively. 

Errors due to relief displacement can be reduced by 
the following methods: 

a. By choosing measurement points that have small 
differences in elevation from the average elevation of the 
terrain. 



b. By choosing measurement points that are 
to the center of the photograph as possible; and 

c. By computing the relief displacement 
measurement point and correcting the photographic 
nates of this point for the amount of the 
displacement . 



as near 

of each 
coordi- 
relief 
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d. 



By making the measurements on a stereo model 
rather on a single photo. 

To calculate the value of the relief displacement, 
the elevation of each measurement point should be found as 
well as the flight altitude, H, and the radial distance, r, 
of the photographic image. The above terms may be applied 
to Equation 2.5 to find the displacement of the image. If 
the displacement, d, is known, the photographic coordinates 
of the point can be corrected. Additional details 
concerning this method will be examined in the following 
section. 

3 . Scale Variation Due to Tilt 

Although photographs for mapping purposes are 
supposed to be taken vertically, often a small angular tilt 
of the aircraft occurs. The tilt angle is usually less than 
two degrees (Moffitt , 1980). The exact scale of a tilted 
photograph depends on the camera focal length, f , airplane 
altitude, H, above datum and the amount of the tilt. Methods 
exist to eliminate scale variations due to tilt, but require 
special instruments. Scale variation due to tilt is not 
eliminated in the present study. All photographs were 
assumed first to be truly vertical. If during the 
measurements the scale variation for a stereo model exceeds 
five percent, then the model is rejected. Five percent scale 
variation gives an error in the measurements of 0.15 ft 
(0.045 m) for each 3 ft (0.91 m) of ground distance. 
Additional details about tilted photographs are discussed in 
Appendix A. 
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C. DESCRIPTION OF RESEARCH PROCEDURE 



1. Procedure to Select , Locate and Purchase Aerial 

Photographs 

The first step was to select a location for study in 
order to select the aerial photographs. Several requirements 
were used to select the most suitable aerial photographs. 
These requirements were: resolution, height of the flight, 
sharpness, accuracy, and maximum scale. Once a photograph 
was selected, the next step was to locate the needed photo- 
grammetric information, i.e., focal length, flying height, 
et c . 

The areas studied were the shoreline and adjacent 
cliff between Monterey Harbor ( 36^ 36' 05"N x 121* 53' 
04"W) and Fort Ord's Stilwell Hall (36* 40’ 00”N x 121* 49' 
06"W) , four specific locations were chosen for the study. 
The chosen areas were significant, because the success of 
the study of coastal erosion was contingent upon the avial- 
ability of suitable aerial photography. 

To select the proper aerial photographs an approxi- 
mate idea of the rate of erosion to measure was needed. 
Several studies of the region (Moffitt, 1968; Thompson, 
1981) indicated average rates of erosion of 3.4 ft/yr (1.04 
m/yr) . On the basis of these it was established that 
photographic resolution for erosion measurements in a 
period of one year should be at least 3 ft (0.91m). 

Another important factor was horizontal "resolu- 
tion", the ability to separate adjacent features so that 
they can been seen as individual images and expressed as 
the maximum number of lines/mm that can be resolved. The 
standard emulsions used in photo gramme try yield a final 
image resolution of 40 lines/mm or more. When copies are 
made from original exposures, each copy of the image can 
degrade the resolution by about 15 to 20 percent. This 
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indicates that the resolution would be 36 lines/mm or more 
for the first copy. For the purpose of this study copies 
were used because the owners of the aerial photographs would 
not allow work with the originals. Therefore, the 
36-lines/mm resolution allowed the use of copies of aerial 
photographs to scales of 1:36000 (Keller, 1975). Sharpness 
of imagery was a prime concern in the choice of the photo- 
graphs. Films such as aerial Ektachrome provide more infor- 
mation to the viewer due to their higher spectral content. 
Black and white has been used for years with success, but 
resolution is limited by silver halide grain size. 

Once criteria had been established for the selection 
of aerial photographs, an attempt was made to obtain 
existing photography, and to choose photographs best suited 
for a coastal erosion survey. Various government agencies 
that maintain permanent records of local aerial photographs 
were contacted. Sources most important in this respect are: 
U.S. Army Corps of Engineers; National Cartographic 
Information Center (NCIC) ; and the University of California 
at Santa Cruz. 

NCIC has the most extensive records to be found. It 
has a system which is called the Aerial Photography Summary 
Record System (APSRS), this system offers a simple method to 
determine not only whether aerial coverage is available over 
a particular geographic location, but who was responsable 
for the photographic project, and the photographic parame- 
ters, i.e. the focal length, photographic scale, and film 
emulsion, etc., associated with each particular roll of 
film. A single set of aerial photographs for 1966 was 
purchased from NCIC. 

The U.S. Army Corps of Engineers, San Francisco, 
indicated that they had sets of aerial photographs from 
1946. We had copies of these printed at the U.S. Geological 
Survey (USGS), Palo Alto. The Corps of Engineers offices 
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in Los Angeles sent a set of copied aerial photographs for 
1984; that flight was made expressly for the purpose of 
studying coastal erosion. 

The Library of the University of California, Santa 
Cruz, has significant holdings of aerial photographs of 
Santa Cruz and adjacent areas; these include aerial photo- 
graphs for the years 1956, 1976,1978 and 1981. The photo- 
graphs are available as paper prints and are easy to locate 
by means of an index. 

For future studies of the area, all the photographs 
we collected are available as strips that cover the entire 
shoreline from Monterey harbor to Fort Ord. Photographs 
purchased from NCIC, were printed with particular attention 
paid to enhancing the cliff tops. The size of the photo- 
graphs used was 9" x 9"; only three sets were available in 
color, while the remaining three sets were in black and 
white . 

2 . Selection and Description of Sub j ect Areas 

The locations of the four areas chosen for study are 
indicated in Figure 2.4 . The general criteria for choosing 
the study areas were that a cliff be adjacent to the shore- 
line, the location be close to a development which had been 
affected by erosion, and that the area had been studied 
previously so that earlier findings could be compared with 
the proposed model. The distribution of the chosen areas 
gives a clear representation of the erosion rates between 
Monterey Harbor and to just north of Fort Ord. 

The area adjacent to the Naval Postgraduate School's 
Beach Laboratory is approximately 1800 ft (549 m) long, and 
the average cliff elevation is approximately 19 ft (5.8 m) 
above Mean Lower Low Water (MLLW) The developments that 
bounded the area are the Monterey Regional Water Pollution 
Control Agency building to the north and the Del Monte Beach 
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Figure 2.4 Selected Areas of Study 

Town Houses to the south. No erosion studies had been 
previously conducted in this area, although the geographic 
situation was ideal for the study. 

The waterfront of the Phillips Petroleum Property is 
approximately 660 ft (201 m) long, and the average cliff 
elevation is approximately 24 ft (7.4 m) above MLLW. The 
development that bounds the northern edge of the area is the 
Holiday Inn, where a vertical concrete wall was built to 
provide protection, and the development that bounds the 
southern edge is the Ocean House Apartments, where a granite 
rock revetment wall was constructed for erosion protection 
in 1983. Studies conducted by Thompson (1981) indicate an 
erosion rate here of 1.8 ft/yr (0.54 m/yr) between 1939 and 
1978 . 

The waterfront of the Monterey Sand Company in Sand 
City is approximately 2200 ft (671 m) long and does not have 
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a development close by. However, it is possible that this 
sand mine contributes to the erosion of the southern area of 
the Monterey Bay (Moff itt , 1968). 

The Fort Ord waterfront in our area of interest is 
approximately 3280 ft (1000 m) long, and the average cliff 
elevation is approximately 52 ft (15.85 m) above MLLW. 
Stilwell Hall is located at the center of this area and is 
in danger due to the erosion of the adjacent cliff. Smith 
(1983) summarized the erosion here as 6-7 ft/yr. 

3 . Selection of Reference Points 

After selecting the four study areas, the next step 
was to choose stable points to reference the location of the 
shoreline. To achieve the greatest accuracy, special atten- 
tion was given to the selection of reference points. 
Objects with sharp, well defined images on the photographs 
were selected as close as possible to the area being studied 
to minimize errors due to scale variation or tilt. The 
height of the reference object is not a problem, because the 
error due to relief displacement can be corrected. For 
example, in the area corresponding to the Beach Lab (NPS), 
street corners between Second Street and Sloat Avenue were 
choosen as reference points. The distances between these 
points were surveyed in the field to determine the scale of 
the photo. 

The geometrical distribution of the reference point 
is important, because to compare identical areas for 
different years, it is necessary to specifically match 
geometrical figures. 

4. Selection of Measurement Points 

Several measurement points along the beach profile, 
such as high tide line, toe of the cliff, and top of the 
cliff, were considered (see Figure 2.5). Past measurements 
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indicate the beaches of southern Monterey Bay are eroding. 
Therefore, it was desired to choose a point that would be 




sharply defined in the aerial photograph and indicative of 
permanent change (erosion) . 

The beach profile at mean waterline is highly vari- 
able. Variability other than that due to erosion may be 
caused by: (1) Sea level changes due to tides (approxi- 
mately 12-hour periods); (2) Sea level changes due to atmos- 
pheric conditions (hours to days variation); and (3) 
Seasonal variations in beach profile due to onshore- of f shore 
sand transport. Similarly, variations at the high waterline 
may be caused by these factors . The bottom of the cliff can 
be affected by factors (2) and (3). Another factor which 
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can affect the measurements is slumping of cliff material, 
causing displacement of the base of the cliff . 

The change in the location of the top of the cliff 
is a measure of permanent erosion; seventy per cent of the 
coast of Southern Monterey Bay is adjacent to cliffs having 
an average height of 30 ft. The erosion of cliff tops is due 
to slumping caused by erosion at the base of the cliff. 

To use aerial photography several other factors must 
be considered which affect the accuracy of the measurement. 
In panchromatic black and white photographs the definition 
of the water line is not clear at scales smaller than 
1/10000, and errors larger than 3 ft (1 m) are easy to make; 
wind surges have the effect of causing the high tide water 
level to be slightly higher or lower than the level due to 
tide alone, which causes some variations in horizontal loca- 
tions of the high water line. A correction or adjustment of 
the high water line location to account for wind tides is 
not possible because of the lack of local wind data. In 
several black and white aerial photographs the bottoms of 
the cliffs were not clearly defined due to the shadow of the 
cliff. Even if the bottom of the cliff could be measured in 
one year, another year could indicate an apparent shoreline 
change due to the collapse of the cliff wall. 

The sharp representation of the top of the cliff in 
aerial photographs offers a clear and identifiable point of 
reference. Therefore, long-term erosion was inferred by 
measuring the locations of the cliff tops. 

To draw the contour of the cliff, reference points 
were selected as near to the beach as possible. The hori- 
zontal choice of points along the cliff was based on the 
topography and profile of the cliff. For instance, where 
the cliff formed a straight line, only two points were 
necessary, but in those cases where the cliff showed a 
different form, a sufficient number of points were measured 
to define the shape. 



37 



5 . Measurement Techniques 



After all reference points and contours of the cliff 
had been selected on the two sets of aerial photographs 
being compared, the measurement process began. The X-Y 
coordinates of photographic image points were obtained using 
a stereo comparator (Carl Zeiss 1944, Model # 96180). 

A comparator is an instrument having two linear 
scales accurately assembled perpendicular to one another. 
The coordinates are read by comparing these linear scales. 
Precision is improved by the use of verniers, or microme- 
ters, and other devices for reading the scales, and magnif- 
iers to enable better pointing at the photographic images; 
thus, the absolute coordinate of images on a photo may be 
read to 0.01 mm. Measurements can be made separately on 
each photo of a pair, but point transfer is much better if a 
pair of comparators are assembled so that the photographs 
may be viewed stereoscopically . The absolute coordinates 
(X,Y) of each point can be obtained directly from the equip- 
ment, and the difference in elevation between two points can 
be determined directly by measuring the difference of 
parallax between them. 

The parallax of a point measured on a pair of over- 
lapping vertical photographs is equal to the X-coordinate of 
the point measured on the left-hand, X^ , photograph minus 
the X-coordinate of the point measured on the right-hand, X^ 
photograph. In this definition, the X-axis passes through 
the principal point and is parallel with the flight line, 
and the Y-axis passes through the principal point and is 
perpendicular to the flight line. Then, 

dp = - X,. (2.7) 

where dp is the parallax deference. 
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To find the difference in elevation using the 
parallax measurement, the following equation was utilized: 

dh = (H dp) / (dp + b ) (2.8) 

where H is the flight height above datum, b is the distance 
between the principal points of the two photographs (air 
base), and dh is the height difference. 

To set up the photographs on the comparator it was 
necessary to determine the principal point of each photo- 
graph. For this purpose, it was necessary to utilize an 
overlay and to draw over the overlay two lines connecting 
opposing fiducial marks in the photograph. The intersection 
of these lines is the nominal principal point. Once the 
principal point of the photograph was established, it was 
transferred visually to the conjugate point corresponding to 
that principal point in the other photograph. Each photo- 
graph then had two points, the principal point and the 
conjugate point. These two points in each photograph were 
connecting with a straight line to draw the air base . The 
air base was aligned with the X-axis of the equipment to 
eliminate Y parallax. 

Measurements were made as follows: 

a. The first step was to find the X- and 
Y- coordinates of the principal points, and using them in 
Equation 2.9 to find the length of the air base. 

(D.a ) 2 = Uo, ~ X„ 2 ) 2 + (Y o, - Y o 2 ) 2 (2.9) 

where is the distance between the principal points or 
air base (b), Xo, , Y o, represent the coordinates of the 
principal point of the left photo. Xo 1 , Yo t represent the 
coordinates of the principal point of the right photo. 
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b. In the second step, the X- and Y- coordinates and 
X-parallax of the reference points were determined; in each 
set of photographs four reference points, A,B,C, and D, were 
always established. The ground distances AB and CD were 
always known by physical measurement in the field; there- 
fore, by calculating the distances ab and cd from the photo- 
graphs the scale of the photographs could be determined 
using Equation 2.4 . The reference parallax of the photo- 
graph was always the parallax of the reference points. These 
points were always chosen at the same elevation; therefore, 
any parallax difference was due to tilt and not to relief 
displacement. Tilt errors were reduced by averaging the 
scales . 

c. The third step was to measure X- and 
Y- coordinates and X-parallax of the photographic images of 
the ground reference points and the points along the shore- 
line. All of the data were recorded on a data sheet. The 
procedure is explained in Appendix B. 

d. The fourth step was to input all recorded data 
into the computer program (Appendix B) to reduce all of the 
photographic coordinates to ground coordinates at scale 
1:2000 and to correct coordinates for relief displacement. 

e. The fifth step was to plot all reduced coordi- 
nates on a chart at a scale 1:2000. For each area we 
obtained one chart per photograph and then matched all 
together at the reference points. For instance, for the 
Beach Lab area the south corner of the intersection of 2nd 
Street and Sloat Street was used as the reference point. 

f. The sixth step was to find the following parame- 
ters: the average rate of erosion, its standard deviation, 
and the average volume eroded using the average elevations. 
This calculation was performed by a computer program using 
linear interpolation. 
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D. APPLICATION AND ANALYSIS OF COASTAL EROSION DATA 



1 . Presentation of Data 

The cliff top contours for years 1946, 1956, 1966, 
1976, 1978, and 1984 are presented for the four beach 
sections. Shoreline recession and erosion rates are deter- 
mined from the cliff top contours. The cliff top contours 
indicate variability in the recession rates. It is assumed 
that the average erosion rate for each area is more repre- 
sentative than the erosion rates at any one point of the 
study area. Recession rates along the shoreline are not 
uniform. By averaging the recession rates at many locations 
along the shoreline the random measurement errors are 
reduced. The average recession rates for each study area 
are computed from the rates of change at number of points 
along the shoreline and are presented in the following 
sections. The standard deviation is also calculated and 
provides a measure of the variability of the erosion rate 
over the measured shoreline. The average volumetric erosion 
rate is calculated by multiplying the average top recession 
and the average cliff elevation. 

2 . Fort Ord 

Fort Ord area has been divided into two regions , 
North and South of Stilwell Hall (Figure 2.6). The erosion 
directly fronting Stilwell Hall has been retarded by 
construction of a revetment, and the results in this region 
are not representative of the typical erosion. Therefore, 
the middle region was not taken into account in the calcula- 
tion of erosion rates . 

The variations of the shorelines for the period 1946 
- 1984 are illustrated in Figures 2.7 and 2.8 in which it 
can be seen that erosion north of Stilwell Hall is more 
intensive than to the south of it. This is because the 
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Monterey Bay 



N 




STILVELL HALL 
(fort Ord) 



Figure 2.6 Fort Ord (Southern Monterey Bay) 

beach to the south appears to be somewhat protected from 
storm waves from the north and west due to the presence of 
the Stilwell Hall revetment. The recession rates, standard 
deviations, and the average eroded volume are shown in 
Tables 1 and 2 . Data are shown for six time increments to 
give a record of change from 1946 until 1984, a total of 38 
years. Of the five time intervals, the 1978-1984 interval 
had the most severe erosion for the North region, and the 
1956-1966 interval for the Southern region. In the analysis 
of erosion during these intervals, the number and direction 
of storms is important. Storms during 1946-1984 period are 
listed in Appendix D. 
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Figure 2.8 Fort Ord (South) 



TABLE 1 

Erosion Rate Results for Fort Ord (North Region) 
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3 . Sand Dune 



The Sand Dune region is shown in Figure 2.9 . The 
variations of the shoreline are illustrated in Figures 2.10 
and 2.11 . The recession rates, standard deviations, and 
the average eroded volumes are shown in Tables 3 and 4 . 
The most severe erosion period is between 1976-1978 (8.8 
ft/yr or 2.7 m/yr) for the southern region and between 
1946-1966 (9.2 ft/yr or 2.8 m/yr) for the northern region. 
The annual rate of erosion during the measured period 
(1946-1984) is 6.3 ft/yr (1.93 m/yr) for southern region and 
6.4 ft/yr (1.96 m/yr) for northern region, which are essen- 
tially the same rates. 




Figure 2.9 Sand Dune (Southern Monterey Bay) 
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TABLE 3 

Erosion Rate Results for Sand Dune (North Region) 
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TABLE 4 

Erosion Rate Results for Sand Dune (South Region) 
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4.. Phillips Petroleum Property 



The Phillips Petroleum property is shown in Figure 
2.12 . Erosion in the southern part of this region is more 
intensive than in the northern part (Figure 2.13). The 
recession rates, standard deviations, and the average 
eroded volumes are shown in the Table 5 . The most severe 
erosion is seen to have occurred in the 1976-1978 interval. 




Figure 2.12 Phillips Petroleum Property 
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Figure 2.13 Phillips Petroleum Property 
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5 . Beach Lab 



The Beach Lab region is shown in Figure 2.14 . The 
variations of the shorelines are illustrated in Figures 2.15 
and 2.16 . The recession rates, standard deviations, and 
average eroded volumes are shown in Tables -6 and 7 . The 
highest rate of erosion occured during the 1976-1978 
interval, both at the northern end (4.2 ft/yr or 1.3 m/yr) 
and at the southern end (3.1 ft/yr or 0.95 m/yr). The 
average annual recession rate between 1946 and 1984 is 1.9 
ft/yr (0.59 m/yr) for the northern and the southern region. 
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Figure 2.15 Beach Lab (North) 
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6 . Comparison with Earlier Studies and Comments 



The results of this and previous erosion studies are 
summarized in Figure 2.17 for Southern Monterey Bay. The 
average erosion rate is a minimum in the southern part and 
increases northward, with a maximum average erosion rate of 
7.3 ft/yr (2.2 m/yr) at the northern Fort Ord region. No 
significant discrepancies were found between the results of 
the present study and those of earlier studies. For the 
southern part of the Beach Lab region Jones (1983) found the 
recession rate to be 1.1 ft/yr (0.33 m/yr). The difference 
between recession rates of the present study and Jones' is 
+0.8 ft/yr (0.24 m/yr). The recession of the cliff toe at 
the Phillips Property was previously studied by Thompson 
(1981); he found a recession rate of 1.8 ft/yr (0.54 m/yr), 
while we found 2.8 ft/yr. An earlier study of the Sand Dune 
region was performed by Moffitt (1968), who observed a 
recession rate of 6.9 ft/yr (2.1 m/yr). The difference 
between the erosion rate he found and ours is 0.6 ft/yr 
(0.18 m/yr). The Fort Ord region was studied previously by 
Jones (1983). The erosion rates he found do not appear to 
be consistent with ours, but they do show - as we do - that 
the northern part is eroded more rapidly than the southern 
part . 

With only a few exceptions the cliff tops always 
recessed shoreward between photographs. This is contrary to 
the findings of earlier investigators (Thompson, Moffitt, 
Jones). The earlier studies indicate in some time inter- 
vals, accretion occurred along sections of the shoreline. 
This is not surprising, because they measured the position 
of the toe of the cliff or of the waterline, which we have 
shown to be unreliable. Thompson measured the toe of the 
cliff, which would indicate accretion in the case when the 
cliff slumps down. He did not remove errors caused by 
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relief displacement and he did not use an accurate photo- 
grammetric instrument (mirror stereoscope) for his measure- 
ments. Moffitt measured the variation of the shoreline 
using as reference the position of the waterline. He 
corrected the location of the waterline for the tidal fluc- 
tuation, but he did not take into consideration the seasonal 
variation and the meteorological effects in sea level varia- 
tion. Jones measured some simple points (one to three 
points) for each region using a mirror stereoscope without 
removing the effect of the relief. 
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SHORELINE LENGTH (FT) 




III. TIDES 



The range of astronomical tides and sea level varia- 
tions caused by oceanic and atmospheric factors must be 
taken into account when studying coastal processes. A major 
factor responsible for episodic erosion events is the simul- 
taneous occurrence of high tides and storm waves. As the 
sea rises and falls with the tide, different elevations of 
the beach become subject to the action of the sea. On a 
tideless coast, the area of beach coming under any partic- 
ular part of the wave action is small in dimension and is 
limited by the size of the waves themselves; the position of 
the break point remains essentially constant for constant 
wave height. On the other hand, if the tidal range is 
considerable, the break point of the waves moves in and out 
with the tides. The effect of swash can be exerted over a 
wide stretch of beach (King, 1959). The rise and fall of 
the sea due to tides is regular and predictable. Some addi- 
tional, usually minor, short term factors that may influence 
the coastal sea level are: changes in atmospheric pressure 
over the ocean surface, changes in average density of the 
sea water column due to temperature and salinity changes , 
and wind set-up or set-down against the coast due to storms 
(Bretschneider , 1980). 

An example of a combination of factors which can cause 
extreme sea levels occurred along the coast of California 
during the winter of 1982-1983. This winter, now commonly 
known as the " El Nino Winter ”, was in many respects the 
most severe storm season in several decades along the 
Pasific Coast of North America (Quiroz, 1983). The extreme 
high tides were a result of warmer than normal mixed upper 
ocean temperatures associated with a strong 2-year El Nino 
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current, a storm surge due to low atmospheric pressure and 
persistent on-shore winds. On the average, the water temper- 
ature is coolest in spring, resulting in lower sea levels, 
and warmest in autumn, resulting in higher sea levels. 
Flick and Cayan (1984) studied the extreme sea levels which 
occurred in the San Diego area. They concluded that until 
about July, the observed sea level was close to normal, but 
by fall and early winter it exceeded the 1960-1978 average 
by up to 0.5 ft. This condition persisted through 1983. Sea 
level finally returned to near normal in December 1983. 

Tides can be predicted by harmonic analysis using tidal 
constituents. Harmonic analysis of tides is based upon the 
assumption that the rise and fall of the tide at any 
locality can be expressed mathematically as the sum of a 
series of harmonic terms having certain relations to astro- 
nomical conditions. The general equation for the height, h, 
of the tide at any time, t, is given by 



h = Ho* T'nHn COS ( dnt - ( k n~ (V 0 + u) n ) 



(3.1) 



where: Ho = mean value of the tidal observations; f = node 
factor of constituent; Hy,= mean amplitude of the constituent 
in degrees per solar hour; = epoch of the constituent at 
t=0 for period of observation; (Vo + u) 4 = value of argument 
of constituent at t=0 for period of observation; n = corre- 
sponds to particular constituent being computed (Schureman, 
1971) . 

The Monterey tide station is located on Municipal Wharf 
No. 2 where the water depth is approximately 6.8 meters 
(Figure 1.1). Both a standard float-well tide gauge and a 
Fisher- Porter , bubbler type, Automatic Digital Recorder are 
used to record tidal data continuously on a strip chart and 
on punched tape, respectively. The records are sent to NOS 
for processing and analysis. 
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The Tidal Predictions Branch of NOS, Rockville, 
Maryland, performed a harmonic analysis of the observed 
tidal data for 13 years, isolating 37 harmonic constituents. 
These constituents are estimates of the periodicities and 
amplitudes of the components of the tidal force. The values 
of 20 constituents whose amplitudes were greater than 0.02 
ft are used to predict the tides in the present study. 

The data for each tidal constituent include: mean ampli- 
tude of the constituent H, local epoch of the constituent 
(Kappa), and the modified epoch of the constituent (Kappa 
prime). The data are listed in Appendix C. Additional 
information required for each constituent are the node 
factor, f, speed of the constituent, a, and the local 
value (Vo+u) which are found using Tables and formulas as 
follows : 

a. i Speed of the constituent, a, (Table 2 of Schureman, 
1971) 

b. Node factor, f, (Table 14 of Schureman, 1971) 

c. Greenwich (Vo + u) , (Table 15 of Schureman, 1971) 

Local(Vo + u) = Greenwich(Vo + u) - (k'- k) (3.2) 

The tides at Monterey are mixed, predominantly semi-diurnal, 
and are composed of two low and two high water levels per 
24.8-hour tidal cycle. An example of the predicted tides is 
illustrated in Figure 3.1 . 

Maixner (1973) and Bretschneider (1980) examined the 
differences between predicted and observed tides at 
Monterey. In their work, predicted heights were subtracted 
from hourly observed heights to yield non- astronomic resi- 
duals, or tidal anomalies. The differences have a frequency 
distribution which closely resembles a Gaussian distribution 
(Figure 3.2). Analysis of hourly sea level observations 
over the 13-year period of record show that 94.5% of the 
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observations lie within +/- 05 ft. of the predicted tide 
and 99.9 % of the deviations lie within +/- 1.0 ft. 
Deviations of sea level from predicted values are the result 
of the complex interactions of atmospheric pressure, wind, 
variations in velocity of longshore currents and changes in 
average density of the water column (Brets chneider , 1980). 
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Figure 3.1 Typical Tidal Cycle for Monterey Bay 
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Figure 3.2 Tidal Anomalies 



IV. SEA AND SWELL RUN-UP 



A. WAVE CLIMATOLOGY 

Wave action is the most important factor in the study of 
beach erosion. A wave climatology for Monterey Bay was 
calculated using deep water wave spectra obtained from Fleet 
Numerical Oceanography Center (FNOC). The deep water waves 
are transformed by refraction and shoaling into shallow 
water. The purpose of this study is to obtain wave charac- 
teristics in shallow water which are directly related to 
beach erosion. 

The wave information is based on a 21-year (1964-85) 
climatology generated using the spectral ocean wave model 
(SOWM) for the entire North Pacific Ocean. SOWM is the 
present operational wave spectral model used daily by the 
Fleet Numerical Oceanography Center to predict waves in the 
Northern Hemisphere. The model has been operational since 
December 1974. Prior to this date a less accurate non spec- 
tral model was used. Because of the recognized importance 
of long-term wave statistics, it was decided to reanalyze 
the wind fields prior to 1974 and re-hindcast the waves 
using the spectral model to form the basis for a wave clima- 
tology. Thus, meteorological data on winds and weather over 
the Northern Hemisphere were reanalyzed and used to compute 
improved wind fields over the ocean. These improved wind 
calculations were, in turn, used to compute wave spectra for 
a grid of 1,575 points for every six hours for twenty 
years. For the period since December 1974 the "nowcasts" of 
the wave spectra based on observed winds every six hours 
have been saved for use in expanding the climatology. 
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The spectral ocean wave model can be described as being 
composed of wave generation, propagation of swell, and 
dissipation components. The growth of wave energy for a 
given wind velocity at grid points is accomplished by a 
modified Miles-Phillips technique. Growth is limited by the 
Pierson-Moskowit z fully-developed spectrum for a given wind 
speed. Directions are obtained through an equation derived 
by the Stereo Wave Observation Project. The wave energy 
spectrum at each grid point is represented by a 15-frequency 
by 12-direction matrix. 

Energy for each frequency and direction is then propa- 
gated as dispersive swell into adjoining grid points. An 
Icosahedral-Gnomic map projection is used for the world’s 
oceans, as the adjoining grid points form the correct great 
circle routes of swell propagation. The hindcasts are for 
the Northern Hemisphere oceans only and do not include swell 
that has propagated from the Southern Hemisphere. 

Dissipation is included in the model to account for the 
effect of opposing winds and losses at high frequencies that 
occur as the swell propagates away from a generation area. 
The SOWM has been verified (Lazanoff and Stevenson, 1975; 
Pierson, 1982) and shown to give reasonable estimates of 
wave spectra (Thornton, 1983). 

B. WAVE TRANSFORMATION FROM DEEP TO SHALLOW WATER 

The waves are transformed by refraction and shoaling as 
they propagate from deep to shallow water. The wave 
velocity varies along the crest of a wave as the wave is 
moving at an angle to depth contours. This variation causes 
the wave crest to bend toward alignment with the contours. 
This bending of the wave is known as wave refraction. 
Refraction has significant influence on the wave height and 
distribution of wave energy along a coast. The changes of 
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wave direction of different parts of the wave results in 
convergence or divergence of wave energy. Refraction 
diagrams were calculated using a version of Dobson's (1967) 
linear refraction program. A topographic grid was generated 
using digital bathymetry provided by NOAA. The depths are 
the average of all historical measured depths in a 6 second 
grid given by Latitude and Longitude and referenced to 
mean-higher-high- water (MHHW) . Since there were some grids 
without values, an interpolated, filled and smoothed bathy- 
metry was generated using a DISSPLA subroutine. The 
mercator grid was then projected to a cartesian grid 
oriented north-south with a 200 meters grid spacing. 

Assuming conservation of energy flux and long-crested 
waves, the wave variance density in shallow water (subscript 
h) is related to the deep water variance density (subscript 
o) by 



S (f,9) = r(f,9) K 2 (f) S (f ,9 ) J (4.1) 

n r s o o v ' 

"l5 _ 

X 

- 12 . 

where f is the wave frequency; 9 , the wave direction; Kr, 
the refraction transfer function; Ks , the shoaling transfer 
function and J, the Jacobian of the transformation. The 
matrix multiplication is indicated by Equation 4.2 

The shoaling effects are due to a change in the group 
velocity, Cg , of the waves which to the first order is a 
function of depth and period. The shoaling transfer func- 
tion is given by 
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(4.3) 



k s< f > ■ 

Assuming energy density is conserved between wave rays, 
the refraction coefficient is given by 






(4.4) 



where b (subscript o) is the wave ray separation in deep 
water and b is the wave ray separation water having a depth 
of 4m . 

Wave refraction diagrams were calculated for 5 direc- 
tions in 30-degree increments from 185.5 to 335.5 azimuth 
and 12 periods (7.5 to 25 . 7 s). Examples of refraction 
diagrams for wave periods of 15 s and 25.7 s are shown in 
Figures 4.1 and 4.2 . The calculated refraction coeffi- 
cients for each study area are shown in Tables 8, 9, 10, and 
11 . Wave refraction studies show the four sites we studied 
to be protected from open ocean waves from the south and 
north, and only waves from the WNW to W S W can deliver 
significant energy. The study showed that, due to wave 
refraction, the open ocean waves impinge upon the shore 
within a narrow range of angles (Thornton, 1983). A value 
of Kr greater than one results in an amplification of wave 
energy whereas Kr less than one results in a decrease in 
energy. Examination of the refraction coefficient tables 
show that the coefficents are mostly less than one and in 
general are greater at Fort Ord and decrease to the south. 
Amplification does occur for some wave components from the 
west at Fort Ord; winter storms are often from the west. 
Amplification occurs for some wave components from the WNW 
(305.5) at Sand Dune and Phillips Property, but the compo- 
nents for other directions are significantly reduced due to 
refraction. No amplification of components occurs at the 
Beach Lab. 
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Figure 4.1 Wave Refraction Diagram (Direction: WNW 
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Figure 4.2 Wave Refraction Diagram (Direction: WNW, Period: 25.70 s) 



By applying refraction and shoaling transformation func- 
tions to the two-dimensional deep-water wave spectrum, the 
shallow- water wave spectra are obtained from Equation 4.1 
Waves having periods shorter than 7.5 seconds were not 
considered. Such short period-waves are locally generated 
in deep water, and experience indicates they are not repre- 
sentative of the short -period , local ly- generated waves at 
the coast . 

Several convenient measures can be used to describe the 
waves. Since the spread of angles of the shallow water wave 
spectra is small, it is convenient to collapse the shallow 
water wave specta into one- dimens ional wave spectra as a 
function of frequency only. 

S h (f) = ^ S h (f,9) A6 (4.5) 



The significant wave height (the average of the highest 
one- third of the waves) corresponds closely to the wave 
height seen by an observer and is a commonly used wave 
statistic. Hs is calculated as a function of the total 
energy of the spectrum. 



H-S = 4 (I S h (f) At') 



t/2 



(4.6) 



The modal period is the period of the spectral peak and is 
representative of the average period of the waves. 

The speed of energy propagation is given by 



C (f) 
g 




2 k h i 
sinn 2kh^ 



(4.7) 



where C is the phase speed, and the wave number k is related 
to the frequency by 



(4.8) 
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The wave power is recognized as an important parameter asso- 
ciated with beach processes. The average power per unit 
area is given by: 




where the energy density of the waves is given by 

V f) - P g S h (f) (4.10) 

C. WAVE SET-UP AND RUN-UP 

Wave set-up is the increase in mean sea level in the 
shoreward direction across the surf zone due to the shore- 
ward transport of momentum by waves. The maximum set-up 
occurs at the shoreward limit of set-up as indicated in 
Figure 4.3 . The wave set-up was measured on a mild sloping 

beach by Guza and Thorton (1981) and found proportional to 
the deep water wave height. Holman and Sallanger (1983) 
extended the measurements for steeper beaches and a much 
wider range of wave heights finding similar results but with 
some correlation with tidal stage. For high tide condi- 
tions, the setup is given by 

n = 0.14H J/O (4.11) 

where H 5 j0 is the significant wave height in deep water. 

Wave run-up refers to the rush of water from broken 
waves up the beach, and is measured by the vertical eleva- 
tion reached above the still-water level. Coincident 

arrival with the highest high tide may cause the biggest 
erosion of the shoreline. The magnitude of the vertical 

run-up can be determined by the characteristics of the 

waves, by the refractive effects of the bottom topography 
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and by the configuration of the beach. Guza and Thornton 
(1981) found that in their measurements on mild beaches the 
vertical extent of run-up, Rv , was proportional to the deep 
water wave height. Again, Holman and Sallanger (1983) 
extended the range of measurements and suggest for high tide 
conditions 



R v = 0.88H XO (4.12) 

The set-up and run-up formulations will be incorporated in a 
simple erosion model described next. 




Figure 4.3 Change in Mean Water Level 
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V. PREDICTIVE BEACH EROSION MODEL 



A. MODEL FOR BEACH PROFILE RESPONSE 

A simple model is proposed here based on the hypothesis 
that erosion only occurs infrequently when the water eleva- 
tion, S, exceeds the toe of the cliff elevation T. Each 
time that the total water level, S, is greater than the toe 
elevation, erosion can occur in the cliff (see Figure 5.1). 




Figure 5.1 Characteristics of the response Model 



The cliff recession is assumed to be given by 

R = k(S'/tan^) (5.1) 
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where, R, is the magnitude of the horizontal recession, $ is 
the slope of the beach, and S’ is the given by 

S ’ = S - T (5.2) 



The slope for each region was calculated using the slope 
between the isobath contours of 0 and 2 fathoms obtained 
from the field sheet DA/NPGS-79. The average slopes are 
shown in Table 12 . The average top and toe elevations of 
the cliff were fcund by leveling in the study areas. The 
average elevations are indicated in Table 12 . The propo- 
tionality factor k accounts for differences in material and 
other (unknown) factors. 

The total water elevation S, is calculated by 



s = + H + R v 



(5.3) 



where, ^iob , accounts for the tide height above MLLW , 
predicted by harmonic analysis (see Chapter 3), the wave 
set-up, h, given by equation 4.11 and the vertical extent 
of wave run-up, R y , given by equation 4.12 . The wave 
set-up and run-up formulas are parameterized on the signfi- 
cant deep water wave height, H . A corresponding deep 
water wave height was calculated by taking the shallow water 
wave height calculated in 4 meters back out into deep water 
without accounting for refraction assuming conservation of 
energy flux as follows: 



H 



.0 



~ H S,4M ( 



K s (f p 



(5.4) 



where fp is the peak frequency in the shallow water wave 
spectrum. The maximum water elevation due to tides and 
waves is predicted hourly for the intervals between photo- 
graphs and the predicted recession using equation 5.1 
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summed. The predicted recession rates were obtained by 
application of the model, using a proportionality coeffi- 
cient k = 0.000096 (see Table 4.11). 

The measured and predicted recession rates for each time 
period (1966-1976, 1976-1978, 1978-1983) normalized by the 
values for the entire period (1966-1983) are compared in 
Figure 5.2 . The north and south recession rates for Ft. 
Ord, Sand Dune and Beach Lab areas are averaged to give 
single measured values. The results show that the model 
gives a reasonable indication of the temporal variability of 
erosion, i.e. the model predicts high erosion when the 
erosion was measured high during a particular time interval 
and vice versa 

The average measured and predicted recession rates for 
the entire comparison interval 1966-1983 are compared in 
Figure 5.3 . The model predicts maximum erosion at Ft. Ord 
decreasing south, but predicts less erosion at Sand Dune 
than at Phillips Petroleum. The high measured recession 
rate relative to the predicted rate at the Sand Dune area 
suggests an anomalously high recession rate. The higher than 
expected recession rate may be associated with the sand 
mining of the beach at this location. 
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5.2 Comparison of Normalized Measured 
and Predicted Recession Rates 
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Figure 5 . 3 



Measured vs. Predicted Recession Rates (1966-83 
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VI. CONCLUSION 



The primary objective of this study is to develop and 
evaluate a model for the prediction of coastal erosion rates 
in Southern Monterey Bay. A sequence of six aerial photo- 
graphs taken during the interval 1946-1984 are used to 
measure coastal erosion. It is demonstrated that the aerial 
photographs provide a suitable means to determine coastal 
changes. Previous studies had not used accepted photograin- 
metric techniques to compensate for scale variation due to 
relief, parallax and plane tilt. Careful analysis was 
applied here to minimize errors, which is discussed in 
detail. Scale variations between photographs were estimated 
by deterimining the average scale of each photograph from 
horizontal ground points. Errors due to relief displacement 
were minimized by measuring the X-parallax of each point. 

It has been shown that the top of the cliff offers a 
consistent measurement point for erosion studies because it 
is affected only by storm waves with simultaneous occurrence 
of high tides. It is found that the erosion along the coast 
of Southern Montery Bay is progressive -- there is no accre- 
tion. The average cliff top recession rate is a minimum in 
the southern part of the Bay and increases northward, with a 
maximum average erosion rate of 7.3 ft/yr (2.2 m/yr) at Fort 
Ord. Topographic maps combined with field surveys show that 
the cliff height decreases to the south. 

A simple predictive model is proposed based on the 
hypothesis that the permanent beach erosion in Southern 
Monterey Bay is episodic, occurring infrequently when high 
tides coincide with stormy weather which allows wave action 
to erode the toe of the cliffs. The tides are predicted 
using harmonic analysis. The wave heights are calculated in 
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shallow-water by refracting deep water directional wave 
spectra provided by Fleet Numerical Oceanography Center 
(FNOC) . Erosion occurs when the total water elevation 
exceeds the toe of the cliff elevation. The total water 
elevation is the combination of tides plus set-up and 
run-up, which are proportional to the wave height. The model 
is calibrated against the measured recession rates for the 
years 1966-1983 when the wave data are available. The model 
quantitatively predicts the temporal variability of the 
recession rates occurring between the aerial photograph 
intervals. The model reasonable predicts the spatial varia- 
tion of the recession rates, indicating maximum erosion at 
Fort Ord and decreasing to the south. The variability in 
erosion is attributed to variations in the incident wave 
energy along the shoreline due to wave refraction. 
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APPENDIX A 



SCALE VARIATION DUE TO TILT 

The concept of a tilted photograph is presented here. 
Some commonly used terms for tilted photographs are illus- 
trated in Figure A.l . 




Figure A.l Elements of Tilted Photographs 

The Nadir point, n, is that point at which a vertical line 
through the perspective center, L, of the camera lens 
pierces the plane of the photograph. Tilt, t, is the angle 
between the vertical line and the optical axis, Lo . The 
isocenter, i, is the point on a photograph intersected by 
the bisector of the tilt angle. The isocenter is signifi- 
cant because it is the center of radiation for displacements 
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of images due to tilt. The axis of tilt is a line perpen- 
dicular to the principal line and passing through the 
isocenter. Swing, s, is the angle at the principal point of 
a photograph which is measured clockwise from the positive 
y-axis to the principal line at the nadir point, n,. 

The effect of tilt in an aerial photograph is to cause 
the scale to vary throughout the picture, even if the ground 
is flat and level. If the scale near the center is correct, 
then the scale is smaller on the side that is tilted upward 
and larger on the side that is tilted downward. 




Figure A. 2 Scale Variation on Tilted Photograph 

This is illustrated in Figure A. 2, where dashed lines 
represent a square grid on a flat and level ground as it 
would appear on a truly vertical photograph. The solid 
lines represent the same grid as it would appear on a photo- 
graph that has been tilted toward the upper part of the 
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diagram. The scale of the grid is the same on both the 
truly vertical and the tilted photograph only along the axis 
of tilt (Moffitt, 1980). 

The scale of any image on a tilted photograph is given 
by 



S = (f - y sint) / (H - h) (A.l) 

where t is the tilt angle, y is the distance of the image 
from the isocenter measured in the direction of the tilt, 
positive on the upper side, and where the other terms have 
already been defined. 

From Equation A.l it is obvious that the scale depends 
on the amount of the tilt, the position of the image on the 
photograph relative to the direction of the tilt and f, H 
and h. Also, the magnitude of the scale variation across a 
tilted photograph is directly proportional to the change in 
the distance y of an image measured in the direction of 
the tilt. For tilts of less than approximately three 
degrees, scale variations are relatively small and can be 
ignored. Tilt is most objectionable where photographs are 
being used to determine elevations of points on the ground 
(SI ama , 1980). 

There are several methods to determine the tilt of a 
photograph. To compute the tilt, ground positions of three 
or more image points in the photograph must known. These 
methods are discussed in Slama (1980). The least compli- 
cated method is to rectify aerial photographs . 
Rectification is the process of projecting the negative of a 
photograph, utilizing a special projection printer, from its 
plane onto another plane by translation, rotation, and scale 
change. The negative is tilted about the two axes by 
computed amounts equal to the tilt of the photograph. 
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In the above process the rectified photographs usually 
are enlarged to adjust all the photographs to a common 
scale. After rectification of the photograph, the scale 
will not be further affected by the tilt. 

Most agencies have specific limits on the amount of 
tilt permitted. The Agricultural Stabilization and 
Conservation Service specifies that a resultant tilt 
exceeding four degrees may be cause for rejection of photo- 
graphs (Stafford, 1971). In case rectified enlargements 
are not available, the simplest method which can be utilized 
to reduce the error caused by scale variation due to tilt is 
to compute several scales by the methods referred in the 
previous section and calculate their average. 
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APPENDIX B 



RECORDING OF DATA AND COMPUTER PROGRAMS 

All measurements were recorded on a specially designed data 
sheet that included all data necessary to process the infor- 
mation. It was designed with the same format as the input 
data program. The following data are included on the 
recording sheet (Figure &57) : (a) Ground distance between 
points A and B and between C and D; (b) scale of the chart; 
(c) reference parallax; (d) focal length of the aerial 
camera; (e) year of the photograph and code of the area, 
i.e. Fort Ord 1978 was coded "1978fo"; (f) X- and Y- coordi- 
nates of the principal point of the left photograph; (g) X- 
and Y- coordinates of the principal point of the right 
photograph; (h)Identif ication of the point, X- and Y- coor- 
dinates, and X parallax of the reference points A, B, C and 
D; (i) identification of the points X- and Y- coordinates, 
and their X- parallaxes. 

Upon completion of the measurement and recording of the data 
necessary to make a comparison of a beach location on two 
sets of photographs, data were processed. The calculations 
necessary to reduce the data were highly repetitive and 
ideally suited for the computer. For this purpose, three 
computer programs were written. The first program trans- 
formed the photo coordinates input data to ground coordi- 
nates at a scale of 1/2000. These coordinates were 
corrected for relief displacement using Equation 1.4; the 
second program plotted the reduced coordinates on a chart of 
dimensions 18" x 36"and the third program using linear 
interpolation to calculate the average rate of erosion, its 
standard deviation, the area eroded, and average volume 
eroded using an average elevation. 
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APPENDIX C 



TIDE TABLES 

Values for each one of the tidal constituents, having 
amplitude greater than 0.02 feet, are listed in Table 14 . 
The calculated terms for each one of the costituents are 
listed in Table 15 . 
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TABLE 14 

Tidal Constituents for Monterey, California 



Monterey Municipal Wharf No. 2 
1 January, 1974 - 31 December, 1974 
Prepared by: 

Tides Branch 
National Ocean Survey 



Constituent 


H 


k 


A 


D 

-k’ 




Amplitude 


Epoch 


k' -k 


M2 


1.625 


297.400 


11.9 


- 309 . 3 


S2 


0.425 


295 . 500 


3.8 


- 299 . 3 


N2 


0.366 


272.000 


16.3 


- 288.3 


K1 


1.216 


97 . 800 


1.6 


- 99.4 


01 


0.763 


81.400 


10.4 


- 91.8 




0.069 


279.400 


15.4 


- 295 . 1 


h . 


0.046 


234.400 


20.0 


- 254.4 


(2N2) 


0.046 


248 . 500 


20.6 


- 269 . 1 


( 001 ) 


0.039 


119.600 


- 7.2 


- 112.4 


J* 


0.011 


296 . 500 


8 . 1 


- 304.6 


SI 


0.038 


202.300 


1.9 


- 204.2 


Ml 


0.117 


114.800 


5 . 9 


- 120.7 


J1 


0.071 


107.000 


- 2.8 


- 104.4 


P' 


0.029 


74.300 


14.1 


00 

00 

1 


Q 1 


0.137 


72.900 


14.7 


- 87 . 6 


T2 


0.025 


295.500 


4.1 


- 299.6 


( 2 Q ) 1 


0.020 


65.000 


19 . 1 


- 84.1 


PI 


0 .381 


92.700 


2.2 


- 94.9 


L2 


0.046 


322.800 


7 . 6 


- 330.4 


K2 


0.121 


287 . 700 


3 . 1 


- 290.8 
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TABLE 

Calculated Constituent 


15 

Terms for 


Monterey 


Monterey Municipal 


Wharf No. 


. 2 




1 January, 1974 - 31 


December , 


1974 




Constituent 


Prepared by 
The Authors 
a 


f 




(V+u ) 




Speed 


Node Factor 


Argument 


M2 


28.9841042 


0.984 




173.4 


S2 


30.0000000 


1.000 




-3.8 


N2 


28.4397295 


0.988 




265.7 


K1 


15.0410686 


1.010 




17.6 


01 


13.9430356 


1.015 




151.6 


^2 


28.5125831 


0.988 




256.0 


h 2 


27.9682084 


0.988 




-11.6 


(2N2 ) 


27.8953548 


0 . 988 




-1.9 


(001) 


16.1391017 


1.049 




71.8 




29.4556253 


0.988 




270.7 


SI 


15.0000000 


1.000 




178.1 


Ml 


14.4966939 


1.274 




152.8 


J1 


15 . 5854333 


1.023 




289.4 




13.4715145 


1.015 




234.0 


Ql 


13.3986609 


1.015 




244.0 


T2 


29.9589333 


1.000 




-1.9 


(2Q)1 


12.3542862 


1.015 




336.3 


PI 


14.9589314 


1.000 




347 . 5 


L2 


29 . 5284789 


1.142 




274.8 


K2 


30.0821373 


1.006 




215 . 1 
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APPENDIX D 



HISTORIC OCEAN STORMS 

Historic ocean storms for Monterey Bay are listed in 
Table 16 . They are collected from earlier studies 
performed by Tompson (1978) and Jones (1983). 
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Historic 


TABLE 16 

Ocean Storms (Monterey Bay) 




1 


January 


,1946 - 31 December, 


1983 


Date- 


Year 


Source 


Direction/ 
Type of Storm 


4 


Mar . 


1946 


Jones 


North Winds 


28 


Jan. 


1947 


Jones 


Norterly Gale 


4 


Apr . 


1947 


Jones 


Northerly 


23 


Feb . 


1948 


Jones 


Northwest 


2-3 


Jan. 


1949 


Jones 


Heavy Winds 


27-29 


Oct . 


1950 


Jones 


Northerly Gale 


6 


Mar . 


1951 


Thompson 


Northwest 


2 


Dec . 


1951 


Jones 


Southeastern 


5 


Dec . 


1951 


Thompson 


Northwest 


7 


Dec . 


1952 


Thompson 


Westerly 


23 


Feb. 


1953 


Jones 


Northeast Gale 


13 


Nov . 


1953 


Jones 


Southerly 


1 


Dec . 


1953 


Thompson 


Northeast 


13 


Feb. 


1954 


Thompson 


Westerly 


19 


Dec . 


1955 


Thompson 


Southwest 


23 


Dec . 


1955 


Thompson 


Southwest 


21 


Feb. 


1956 


Thompson 


Westerly 


4 


Mar. 


1956 


Thompson 


Westerly 


3 


Apr . 


1958 


Thompson 


Northwest 


14 


Nov . 


1958 


Thompson 


Northwest 


9-10 


Feb . 


1960 


Jones 


Southerly and 


1 


Feb. 


1963 


Jones 


Westerly 

Northwest 
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Historic 


Ocean 


Table 16 

Storms (Monterey 


Bay) (cont 'd. ) 




7 


Mar . 


1964 


Thompson 


Northwest 




28 


Dec . 


1965 


Thompson 


Southwest 




24 


May . 


1967 


Thompson 


Northwest 




3 


Dec . 


1967 


Thompson 


Northwest 




10 


Dec . 


1968 


Thompson 


Northwest 




6 


Feb . 


1969 


Jones 


- 




24 


Feb. 


1969 


Thompson 


Southwest 




11 


Dec . 


1969 


Thompson 


Southerly 




15 


Jan . 


1971 


Thompson 


Southerly 




22 


Jan . 


1972 


Thompson 


Southwest 




16 


Jan . 


1973 


Thompson 


Southwest 




24 


Feb. 


1973 


Thompson 


Southerly 




28 


Mar . 


1973 


Thompson 


Northwest 




16 


Jan . 


1974 


Thompson 


Southwest 




1 


Mar . 


1974 


Thompson 


Southwest 




9 


Jan . 


1978 


Jones 


South and 




16 


Jan . 


1978 


Jones 


Southwest 

Southwest 


25 


-27 


Jan . 


1978 


Jones 


Southwest 




5 


Mar . 


1978 


Jones 


- 


16 


-19 


Feb. 


1980 


Jones 


Southerly 




1 


Dec . 


1982 


Jones 


Northwest 




2 


Dec . 


1982 


Jones 


- 




23 


Dec . 


1982 


Jones 


Winds up to 




28 


Jan . 


1983 


Jones 


70 mi . /hr 
Southerly 




1 


Mar. 


1983 


Jones 


Southwest 
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